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ABSTRACT 

Context. Starbursts are one of the main contributors to the chemical enrichment of the interstellar medium. However, mechanisms 
governing the interaction between the recent star formation and the surrounding gas are not fully understood. Because of their a priori 
simplicity, the subgroup of H ii galaxies constitute an ideal sample to study these mechanisms. 

Aims. A detailed 2D study of the central region of NGC 5253 has been performed to characterize the stellar and ionized gas structure 
as well as the extinction distribution, physical properties and kinematics of the ionized gas in the central ~210 pcxl30 pc. 
Methods. We utilized optical integral field spectroscopy (IPS) data obtained with FLAMES. 

Results. A detailed extinction map for the ionized gas in NGC 5253 shows that the largest extinction is associated with the prominent 
Giant Hii region. There is an offset of ~0'.'5 between the peak of the optical continuum and the extinction peak in agreement with 
findings in the infrared. We found that stars suffer less extinction than gas by a factor of ~0.33. The [S ii]/16717/[S ii]/t6731 map shows 
an electron density (Ne) gradient declining from the peak of emission in Ha (790 cm"^) outwards, while the argon line ratio traces 
areas with ~ 4200 - 6200 cm"''. The area polluted with extra nitrogen, as deduced from the excess [Nii]/16584/Hq', extends up to 
distances of 3'.'3 (~60 pc) from the maximum pollution, which is offset by ~r.'5 from the peak of continuum emission. Wolf-Rayet 
features are distributed in an irregular pattern over a larger area (~ 100 pcxlOO pc) and associated with young stellar clusters. We 
measured He^ abundances over most of the field of view and values of He^^/H^<0.0005 in localized areas which do not coincide, in 
general, with the areas presenting W-R emission or extra nitrogen. The line profiles are complex. Up to three emission components 
were needed to reproduce them. One of them, associated with the giant Hii region, presents supersonic widths and [Nii]/i6584 and 
[S ii]/t/167 17,6731 emission lines shifted up to 40 km s"' with respect to Ha. Similarly, one of the narrow components presents offsets 
in the [Nii]/16584 line of <20 km s"'. This is the first time that maps with such velocity offsets for a starburst galaxy have been 
presented. The observables in the giant Hii region fit with a scenario where the two super stellar clusters (SSCs) produce an outflow 
that encounters the previously quiescent gas. The south-west part of the FLAMES IFU field is consistent with a more evolved stage 
where the star clusters have already cleared out their local environment. 

Key words. Galaxies: starburst — Galaxies: dwarf — Galaxies: individual, NGC 5253 — Galaxies: ISM — Galaxies: abundances 

— Galaxies: kinematics and dynamics 



1. Introduction 

Starbursts are events characterized by star-formation rates much 
higher than those found in gas-rich normal galaxies. They are 
considered one of the main contributors to the chemical en- 
richment of the interstellar medium (ISM) and can be found in 
galaxies covering a wide range of masses, luminosities, metal- 
licities and interaction stages such as blue compact dwarfs, nu- 
clei of sp iral galaxies, or (Ultra)luminous Infrared Galaxies (see 
IConti et al. 2008, and references therein). 

A particularly interestin g sub s et are the H ii galaxies, iden- 
tified for the first time by Hard (Il956h : gas-rich, metal poor 
(l/40Zo<Z<l/3 Zq) dwarf systems characterized by the pres- 
ence of large ionized Hii regions that dominate their optical 

Send offprint requests to: A. Monreal-Ibero 

* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile (ESO 
Programme 078.B-0043). 



spectra (see iKunth & Ostlinll2000l for a review of these galax- 
ies). These systems are a priori simple, which makes them the 
ideal laboratories to test the interplay between massive star for- 
mation and the ISM. 

NGC 5253, an irregular galaxy l ocated in the C entarus A 
/ M 83 galaxy complex ( Karachent sev et all |2007|) . is a lo- 
cal example of an H ii galaxy. This galaxy is suffering a burst 
of star formation which is be lieved to have been triggered 
by an encounter with M 83 dvan den Berghl Il980h . This is 
supported by the existence of the Hi plume extending along 
the optical minor axis \yhich is best explained as tidal debris 
(iKobulnickv & SkiIImanll2008l) . 

NGC 5253 constitutes an optimal target for the study of the 
starburst phenomenon. On the one hand, its proximity allows a 
linear spatial resolution to be achieved that is good enough to 
study the details of the interplay between the different compo- 
nents (i.e. gas, dust and star clusters) in the central region. On 
the other hand, this system has been observed in practically all 
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Table 1. Basic data for NGC 5253 



Table 2. Observation log 



Parameter 


Value 


Ref. 


Name 


NGC 5253 


(a) 


Other designations 


ESO 445- G004, Haro 10 


(a) 


RA (J2000.0) 


13h39m55.9s 


(a) 


Dec(J2000.0) 


-31d38m24s 


(a) 


z 


0.001358 


(a) 


D{Mpc) 


3.8 


(b) 


scale (pc/") 


18.4 




niB 


10.78 


(c) 


Mb 


-17.13 


(c) 


U-B 


-0.30 


(c) 


B-V 


0.50 


(c) 


V-R 


0.32 


(c) 




1.4 X lO** Mo 


(d) 




~0.3(*> 


(e) 




8.95 


(f) 


Iog(L;,/Lo) 


9.21 


(f) 



Grating 


Spectral range 


Resolution 




Airmass 




(A) 




(s) 




L682.2 


6438-7 184 


13 700 


5 X 1 500 


1.75-1.13 


L479.7 


4501-5 078 


12000 


5 X 1 500 


1.11-1.01 



•> We assumed 12 + log(0/H)o = 8.66 lAsplund et al.ll2004h . 

^> NASA/IPAC Extragalactic Database (NED), 
b) - 

d) 

e) 



Sakaietal. (200^ 
Taylor etaU (.200: 



)4 
0^^ 



Kobulnickv & Skillmaij(l2008h . 
Kobulnickvet 31.71999^ 



On account of their irregular structure, a proper character- 
ization of the physical properties of Hii galaxies, necessary to 
explore the interplay of mechanisms acting between gas and 
stars, requires high quality two-dimensional spectral informa- 
tion able to produce a continuous mapping of the relevant quan- 
tities. Such observations have traditionally been done in the op- 
tical and near-i nfrared by mapping the galaxy under study with 
a long -slit (e.g. IVilchez & Iglesias-Paramolll998l: IWalsh & Rovl 
Il989h . This is, however, expensive in terms of telescope time 
and might be affected by some technical problems such as mis- 
alignment of the slit or changes in the observing conditions 
with time. The advent and popularization of integral field spec- 
troscopy (IFS) facilities, able to record simultaneously the spec- 
tra of an extended continuous field, overcomes these difficul- 
ties. Nevertheless, work based on this technique de voted to the 
study of Hn galaxies is still relatively rare (e.g. La2os et alj 



, ', , , , ,. , ,, stuay or Hn galaxies is still relatively rare (e.g. Lagos et e 

S""^"'"^"^"'- ^2201). Re-scaled to the distance adopted here. | 20ollBordalo et d]l2009Hjames et aLlTOllfehrir et aL„ 200' 



spectral ranges from the X-ray to the radio, and therefore a large 
amount of ancillary information is available. 

The basic characteristics of this galaxy are compiled in 
Table [1] Its stellar content has been widel y studied and mor e 
than 300 stellar clusters have been detected dCresci et al.ll2005l) . 
Multi-band photometry with the WFPC2 has revealed that those 
in its central region present typical masses of ~ 2 - 1 20 X 10^ Mp 
and a re very young, with ages of ~ 1 - 12 Myr (e.g. lHarris et alJ 
|2004 . In particular, HST-NICMOS images have revealed that 
the nucleus of the galaxy is made out of two very massive (~ 
1-2x10'' Mq) super stellar cluste rs (SSCs), with ages of abou t 
~ 3.5 Myr, separated by ~ 0'.'4 dAlonso-Herrero et alJ l2004 . 
and which are coincident with the double radio nebula detected 
at 1.3 cm (iTurner et al.ll200(il) . Also, detection of spectral fea- 
tures characteristic of Wolf-Rayet (W-R) stars i n specific region s 
of the galaxy have been reported (e.g. Schaerer et alj |I997|) . 
Recently, seven supernova remnants have been detected in the 
central re gion of t his galaxy by means of the [Fen]/il.644yum 
emission jLabrie & Pritchet 2006). 

NGC 52 53 presents a filamentary structure in Ha (e.g. 
lMartinl[T998h associated with extended diffuse emission in X- 
ray which can be explained as multiple superbubbles around its 
OBs associations and SSCs that are the results of the combined 
action of stellar winds an d supernovae (Strickland & Stevensi 
[19991 [Summers et alJl200l . 

The measured metallicity of this galaxy is relatively low (see 
Table[Tli and presents a generally uniform distribution. However, 
an increase in the abundance of nitrogen in the central region of 
~ 2 - 3 times the mean has been reported (|W alsh & Rovl ll989t 
iKobulnickv et al.lI997l:lL6pez-Sanchez et alj20 07). No other el- 
emental species appears to present spatial abundance fluctua- 
tions. The reason for this nitrogen enhancement has not been 
fully clarified yet although a connection with the W-R popula- 
tion has been suggested. 



barcfa-Lorenzo et al.ll2008l:llzotov et al.ll2006l)" 

Here, we present IFS observations of the central area of 
NGC 5253 in order to study the mechanisms that govern the 
interaction between the young stars and the surrounding ionized 
gas. The paper is organized as follows: section |2] contains the 
observational and technical details regarding the data reduction 
and derivation of the required observables; section [3] describes 
the stellar and ionized gas structure as well as the extinction dis- 
tribution and the physical and kinematic properties of the ion- 
ized gas; section ID discusses the evolutionary stage of the gas 
surrounding the stellar clusters, focusing on the two most rel- 
evant areas of the field of view (f.o.v.). Section |5] itemises our 
results and conclusions. 



2. Observations, data reduction and line fitting 

2.1. Observations 

Data were obtained with th e Fibre Large Array M ulti Element 
Spectrograph, FLAMES (iPasquini et alJ l2002h at Kueyen, 
Telescope Unit 2 of the 8 m VLT at ESO's observatory on 
Paranal, on February 10, 2007. The central region of the galaxy 
was observed with the ARGUS Integral Field Unit (IFU) which 
has a field of view of 1 1'.'5 x 7'.'3 with a sampling of 0.52"/lens. 
In addition, ARGUS has 15 fibers that can simultaneously ob- 
serve the sky and which were arranged forming a circle around 
the IFU. The precise covered area is shown in Figure [1] which 
contains the FLAMES field of view over-plotted on an HST B, 
Hq', I colour image. 

We utilized two different gratings in order to obtain informa- 
tion for the most important emission lines in the optical spectral 
range. Data were taken under photometric conditions and see- 
ing ranged typically between 0'.'8 and I'.'O. The covered spectral 
range, resolving power, exposure time and airmass for each con- 
figuration are shown in Tabled In addition to the science frames, 
continuum and ThAr arc lamps exposures as well as frames for 
the spectrophotometric standard star CD-32 9927 were obtained. 
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100 pc 



Fig. 1. False colour image in filters F435W (B, blue), F658N 
(Ha, magenta), and F8 14W (/, green) for NGC 5253 using HST- 
ACS images (programme 10608, P.I.:Vacca) with the area cov- 
ered by our FLAMES data over-plotted as a white rectangle. The 
orientation and scale for a distance of 3.8 Mpc are indicated. 



2.2. Data reduction 

The basic reduction steps for the FLAMES data were performed 
with a combination of the pipeline provided by ESO (version 
L0i3 via esorex, version 2.0.2 and some IRAlQ routines. First 
of all we masked a bad column in the raw data using the task 
fixpix within IRAF. Then, each individual frame was pro- 
cessed using the ESO pipeline in order to perform bias subtrac- 
tion, spectral tracing and extraction, wavelength calibration and 
correction of fibre transmission. 

Uncertainties in the relative wavelength calibration were es- 
timated by fitting a Gaussian to three isolated lines in every spec- 
trum of the arc exposure. The standard deviation of the central 
wavelength for a certain line gives an idea of the associated error 
in that spectral range. We were able to determine the centroid of 
the lines with an uncertainty of ~0.005A, which translates into 
velocities of ~0.3 km s~'. The spectral resolution was very uni- 
form over the whole field-of-view with values of 0. 178 +0.004 A 
and 0.241 + 0.009 A, FWHM for the blue and red configuration 
respectively, which translates into o-i„stru ~ 4.7 km s"'. 

For the sky subtraction, we created a good signal-to-noise 
(S/N) spectrum by averaging the spectra of the sky fibres in 
each individual frame. This sky spectrum was subsequently sub- 
tracted from every spectrum. In several of the sky fibres, the 
strongest emission lines, namely [Oiii]/I5007 in the blue frames 
and Ha in the red frames, could clearly be detected. We at- 
tributed this effect to some cross-talk from the adjacent fibers. 
A direct comparison of the flux in the sky and adjacent fibers 
showed that this contribution was always <0.6% which is neg- 



' http://www.eso.org/projects/dfs/dfs-shared/web/vlt/vlt-instrument- 
pipelines.html. 

^ The Image Reduction and Analysis Facility IRAF is distributed by 
the National Optical Astronomy Observatories which is operated by the 
association of Universities for Research in Astronomy, Inc. under coop- 
erative agreement with the National Science Foundation. 



ligible in terms of sky subtraction. However, in order to reduce 
this contamination to a minimum, we decided not to use these 
fibres in the creation of the high S/N sky spectra. 

Regarding the flux calibration, a spectrum for the calibra- 
tion star was created by co-adding all the fibers of the standard 
star frames. Then, a sensitivity function was determined with the 
IRAF tasks standard and sensfunc and science frames were 
calibrated with calibrate. Afterwards, frames corresponding 
to each configuration were combined and cosmic rays rejected 
with the task imcombine. As a last step, the data were refor- 
matted into two easier-to-use data cubes, with two spatial and 
one spectral dimension, using the known position of the lenses 
within the array. 

2.3. Line fitting and map creation 

In order to obtain the relevant emission line information, line 
profiles were fitted using Gaussian functions. This procedure 
was done in a semi-automatic way using the IDL based routine 
MPFITEXPrQ (Markwardt 2009) which offers ample flexibility 
in case constraints on the parameters of the fit are included, such 
as lines in fixed ratio. The procedure was as follows. 

As a first step, we fit all the lines by a single Gaussian. The 
Hot -I- [Nil] complex was fitted simultaneously by one Gaussian 
per emission line plus a flat continuum first-degree polynomial 
using a common width for the three lines and fixing the sepa- 
ration in wavelength between the lines according to the redshift 
provided at NEEq and the nitrogen line ratio (/i6583//i6548) to 
3. The same procedure was repeated for the [S ii]/l/l6717,6730 
doublet, the [Ariv]/14711 line (which was fitted jointly with the 
[Feiii]/14701 and He 1/14713) and the [Ariv]/}4740 line (which 
was fitted jointly with the [Feiii]/14734 line), but this time with- 
out any restriction on the line ratios. Finally, H/3, [Oiii]/I5007, 
and He 1/16678 were individually fitted. 

This single Gaussian fit gave a good measurement for the 
line fluxes and the results of these fits are the used in all the 
forthcoming analysis, with the exception of the kinematics. This 
latter analysis requires a more complex line fitting scheme, since 
several lines showed signs of asymmetries and/or multiple com- 
ponents in their profiles for a large number of spaxels. In those 
cases, multi-component fits were performed. Over the whole 
field of view we compared the measured flux from performing 
the fit with a single Gaussian to the fit by several components 
in the brightest emission lines (namely: [Oiii]/I5007, Ho-, 
[N ii]/I6584 and [S ii]/I/I6717,673 1). Differences between the two 
sets of line fluxes ranged typically from 0% to 15%, depending 
on the spaxel and the emission line, and translated into differ- 
ences in the line ratios <0.06 dex. 

In all the cases, MPFITEXPR estimated an eiTor for the fit 
using the standard deviation of the adjacent continuum. Those 
fits with a ratio between line flux and error less than three were 
automatically rejected. The remaining spectra were visually in- 
spected and classified as good or bad fits. 

Finally, for each of the observables, we used the derived 
quantity together with the position within the data-cube for each 
spaxel to create an image suitable to be manipulated with stan- 
dard astronomical software. Hereafter, we will use both terms, 
map and image, when referring to these. 



^ See http://purl.com/net/mpfit 

* See ht tp://nedwww.ipac. caltech.edu/l 
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Table 3. Main reference clusters. 



Continuum 



Name 


FLAMES 




Other names 






coord.(",") 


C97(^) 


H04<'" 




AH04"'> 


#1 


(3.6,0.5) 


N5253-5 


1 


UV3 


C1+C2 


#2 


(LO,-LO) 


N5253-4 


4,8,24,25 


UVl 




#3 


(-4.2, LO) 


N5253-3 


3,5 




C4+C5 



Calzettiet ari(ll997h . 
Harris et alT(l2004h. 



Kobulnickv et alj [l997l) 



Alonso-Herrero et alj ( |2004|) 



3. Results 

3. 1 . Stellar and ionized gas structure 

Figure |2] displays the stellar structure, as traced by a continuum 
close to Ha, as well as the one for the ionized gas (traced by the 
Ha emission line). The over-plotted contours, which represent 
the HST-ACS images in the F659N and F8 14W bands convolved 
with a Gaussian to match the seeing at Paranal, show good cor- 
respondence between the images created from the IFS data and 
the HST images (although obviously with poorer resolution for 
the ground-based FLAMES data). A direct comparison of these 
maps shows how the stellar and ionized gas structure differs. 

The continuum image displays three main peaks of emission 
which will be used through the paper as reference. We have asso- 
ciated each of these peaks with one or more star clusters by direct 
comparison with ACS images. Table [3] compiles their positions 
within the FLAMES field of view together with the names of 
the corresponding clusters according to several reference works. 
These clusters trace a sequence in age as we move towards the 
right (south-west) in the FLAMES field of view. The clusters as- 
sociated with peak #1 are very yo ung (~ 3 - 8 Myr, iHarris et alJ 
l2004t lAl onso-Herrero et al] l2004i) . those associated with peak 
#2 display a ran ge of ages from ve ry young to intermediate age 
(~ 6 - 170 Myr, Harris et al. 2004), while the stars in the pair of 
clusters associ ated with pe ak #3 seem to have intermediate ages 
(70-1 13 Myr, iHarris et al.l l2004). 

The H a emission line reproduces the structure described by 
ICalzetti et al. ( 1997) using an HST WFPC2 image. Briefly, the 
central region of NGC5253 is divided into two parts by a dust 
lane that crosses the galaxy along the east-west direction (from 
~[2'.'0,-3'.'0f| to the Complex #3). Most of flie Ha emission is 
located towards the north of this lane where there is giant H ii 
region associated with the Complex #1. This region shows two 
tongue-shaped extensions towards the upper and lower part of 
the FLAMES field of view (P.A. on the sky of ~ 45° and ~ 
-135°, respectively) as well as a extension at P.A.~ 155° which 
contains the Complex #2. Towards the south of the dust lane, the 
emission is dominated by a peak at ~[-3"5,-2"5 1 which could be 
associated with cluster 17 in IHarris et al.l (l2004l) . 

3.2. Extinction structure 

Extinction was derived assu ming an intrinsic B almer emission 
fine ratio of Ha/Hy8 = 2.87 (lOsterbrock & Ferla nd 2006 , for an 
Te - 10000 K) and using the extinction curve of .Fluks et al.l 
Since the Ha and Hy6 emission lines are separated by 
~ 1700 A and both of them had high signal-to-noise ratio in only 
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^ Hereafter, the different quoted positions will be refereed as [" 
and using the FLAMES f.o.v. as reference. 
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Fig. 2. Top: Stellar component distribution as traced by a con- 
tinuum map made from the average flux in the spectral range 
6525 - 6545 A and 6600 - 6620 A. Bottom: Ionized gas distribu- 
tion as traced by the Ha emission line. We have over-plotted con- 
tours corresponding to the HST-ACS images in the F814W (top) 
and F658N (bottom) filters (programme 10609, PL: Vacca) con- 
volved with a Gaussian of 0'.'8 to simulate the seeing at Paranal. 
The position of the three main peaks of continuum emission are 
marked with crosses. Both maps are presented in logarithmic 
scale in order to emphasize the relevant morphological features 
and cover a range of 2.0 dex for the Ha and 0.9 dex for the con- 
tinuum map. Flux units are arbitrary. 



a single exposure, we decided to obtain the extinction maps from 
the LR3 and LR6 exposures observed at the smallest airmass 
(1.009 and 1.160, respectively), thus minimizing any effect due 
to differential atmospheric refraction. 

We have not included any correction for an underlying stellar 
population. We inspected carefully each individual spectrum to 
look for the presence of the stellar absorption feature in H/3. Only 
in those spaxels associated with the area around the Complex #3 
(~9-10 spaxels or ~l'.'5xl'.'5 in total) was such an absorption de- 
tected. We estimated the influence of this component by fitting 
both the absorption and the emission component in the most af- 
fected spaxel. Equivalent widths were ~3 and ~9 A, respectively. 
In these spaxels, the absorption line was typically ~ 15-20 times 
wider and with about half - one third of the flux of the emission 
line. This implies an underestimation of HyS emission line flux 
of about 10%. For the particular area around Complex #3, this 
translates into a real extinction of ~ 0.9 mag instead of the 
measured Ay ~ 1.2 mag. 
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Fig. 3. Reddening map obtained from the hydrogen recombi- 
nation lines assuming Ha/H/3 - 2.87 and th e extinction law of 
Fluks et al.1 (0994) and E{B - V) = Ay/3.1 jRieke & Lebofskvl 
1985h . The position of the three main peaks of continuum emis- 
sion are shown for reference. We have over-plotted the F160W 
band NICMOS-HST image from lAlonso-Herrero et alj (|2002|) 
with contours for comparison. 



The correspondin g reddening map was dete rmined assuming 
E{B - y) = Av/3.1 fRiek e & Lebofskvlll985l) and is presented 
in Figure [3] This map shows that the extinction is distributed in 
a non-uniform manner ranging from E(B - V) = 0.16 to 0.64 
(mean 0.33, standard deviation, 0.07). Given that G alactic red- 
dening for NGC 5253 is 0.056 ( Schlegel et aO 19981) . nearly all 
of the extinction can be considered intrinsic to the galaxy. 

In general terms, the structure presented in t his ma p co- 
incides with the one presented by ICalzetti et alJ (Il997h . The 
dust lane mentioned in the previous section is clearly visi- 
ble here and it causes extinction of Ay ~ 0.10 - 1.13 mag. 
However, the larger measured extinction values are associated 
with the giant H 11 region, in agree ment with the H i distribu- 
tion dKobulnickv & Skilknanl 120081) . Dust in this area forms an 
S-shaped distribution with Ay ~ 0.13 - 1.15 mag in the arms. 

In order to explore the relation between the extinction suf- 
fered by the gas and by the stellar populations, our E{B - V) 
measurements were compared with colours defined ad hoc. For 
the covered spectral range, it is not possible to exactly simulate 
any of the existing standard filters. It is possible, however, to 
create filters relatively similar to the g' and Rc ones. We have 
simulated two set of filters. 

In the first case, the flux was integrated over two large wave- 
length ranges (465 - 495 nm and 643 - 673 nm) in order to sim- 
ulate broad filters. The relation between the reddening derived 
for the ionized gas and the derived colour, hereafter (^j - Rd), 
is shown in the upper panel of Figure |4l as would be observed 
with photometry. The first order polynomial fit to the data and 
the Pearson correlation coefficient are included on the plot. Also 
shown is the expected relation for an Im galaxy with foreground 
reddening. This latter was derived for the average of two Im tem- 
plates (NGC 4449 and NGC 4485) from Kennicutt (1992) ap- 
plying a fore ground screen of du st with a standard Galactic red- 
dening law ( Cardelli et a n il989h with R=3.1. There is a strong 
difference between the expected relation for a foreground screen 
of dust and the measured values. On the one hand, colours are 
much redder On the other hand, the slope of the 1 -degree poly- 
nomial fit is much less steep than the expected one. Also, there is 
a very good correlation between the E{B - V) and our synthetic 
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Fig. 4. Relation between the derived reddening and a colour sim- 
ilar to g' - Rc- The vertical line represent the g' - Rc colour 
expected for a Im type galaxy without any extinction while the 
green dashed line corresponds to the expected relation if gas and 
stars were suffering the same amount of extinction (see text for 
details). The 1-degree polynomial fit to all the data appear as 
a continuous line. The corresponding fit is shown in the upper 
left corner while the Pearson's coefficient of correlation is indi- 
cated in the lower right corner The simulated g' - Rc colours 
have been done by integrating the flux in the spectral ranges of 
4 653 - 4 953 A (g'j) and 6 43 1 -6 73 1 A (/?,i ) for the upper panel 

and 4 600 - 4 660 A (g',) and 6 600 - 6 680 A (Rci) for the lower 
one. 
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Fig. 5. Line free colour map. The simulated filters have been de- 
fined as explained in Figure |4] The position of the three main 
peaks of continuum emission are shown for reference. 



(g' - Rci) colour. All this can be attributed to the contamination 
of the gas emission lines, mainly Ha and Hy6, in our filters. 

In the second set, we restricted the spectral ranges for the 
simulated filters to a narrower wavelength range which was free 
from the contamination of the main emission lines. The map 
for this line-free colour is displayed in Figure |5] The structure 
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resembles the one presented in Figure [3] (i.e. dust lane, redder 
colours associated with the giant H ii region), although there are 
differences, that can be attributed to differences in the properties 
of the stellar populations in the different clusters. The relation 
between the reddening and the corresponding ig'2-Rc2) is shown 
in the lower panel of Figure |4] This time colours are more simi- 
lar to what is expected for a given stellar population suffering a 
certain amount of extinction. However, for a given colour, stars 
do not reach the expected reddening if gas and stars were suf- 
fering the same extinction (i.e. data points are below the green 
line). The ratio between the slopes indicates that extinction in 
the stars is a factor 0.33 lower than the one f or the ionized gas. 
This is similar to what ICalzet ti et al.l (1 19971) found using HST 
images who estimated that the extinction suffered by the stars is 
a factor 0.5 lower than for the ionized gas and can be explained 
if the dust h as a larger covering factor for the ionized gas than 
for the stars (ICalzetti et al .1119941) . 

In general, our E{B - V) measurements agree with pre- 
vious oiies_jisiiig_flie_sameM areas 
(e.g. iGonzalez-Riestra et al.1 Il987t iL opez-Sanchez et al.l l2007h 
or with poorer spatial resolution Imish & Roy 1989). However 
there are discrepancies when comparing with the estimation of 
the extinction at other wavelengths. In particular, the peak of 
extinction (Ay - 2.1 mag, according to the Balmer line ra- 
tio) is offset by ~ U. 'S from the peak of continuum emission. 
lAlonso-Herre ro et al.l ^2004) showed how in the central area 
of NGC 5253 there are two massive star clusters, CI and C2. 
While CI is the dominant source in the optical, coincident with 
our peak in the continuum map, the more massive and extin- 
guished C2 is the dominant source in the infrared. The con- 
tours for the NICMOS F160W image in Figure [3] show the 
good correspondence between our maximum of extinction and 
C2. Measurements in the near and mid-infrared sugge st ex- 
tinctions of Ay ~ 17 ma^ for this cluster (Tu rner et al.l 20031: 
lAlonso-Herrero et al.l l2004: Mar tin-Hernandez et al. 2005). The 
discrepancy between these two values indicates that a fore- 
ground screen model is not the appropriate one to explain the 
distribution of the dust in the giant H ii region. 



3.3. Electron density distribution 

Electron density (Ne) can be determined from the 
[Sii]^6717/[Sii]^6731 and [Ariv]^4711/[Ariv]^4740 line 
ratios; values of 1.25 and 1.00 respectively were measured in 
a spectrum made from the sum of all our spaxels. Hereafter, 
we will refer to this spectrum as the integrated spectrum. 
Electron densities were determined assuming an electron 
temperature of 1 1 650 K, t he average of the values given in 
iLopez-Sanchez et a n (I2007h. and using the task temden, based 
on the fivel program (IShaw&Dufou3ll995h included in the 
IRAF package nebular. Derived values of A^^ for the two line 
ratios were 180 cm and 4520 cm respectively. Differences 
between the electron densities derived from the argon and 
sulphur lines are usually found in ionized gaseous nebulae (see 
IWang et ani2004l) and are understood in terms of the ionization 
structure of the nebulae under study: [Ariv] lines normally 
come from inner regions of higher ionization degree than [S n] 
lines. Typically, for giant Galactic and extragalactic H ii regions, 
derived Np from the s e two hne ratios diff er in a factor of ~5 (e.g. 
lEsteban et al.l l2002t iTsamis et alj l2003h which is much lower 
than what we find for the integrated spectrum of NGC 5253 
(~25). However, when only the giant Hii region is taken into 
account (i.e. the area of ~90 spaxels where the argon lines are 
detected) the difference between the densities derived from the 
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Fig. 6. Maps for the line ratios sensitive to the electron density. 
The position of the three main peaks of continuum emission 
are shown as crosses for reference. The displayed ranges in the 
line ratios imply electron densities of < 10-790 cm"^ and 190- 
8750 cm"^ for the sulphur and argon line ratio, respectively. 



argon and the sulphur lines (~10, see typical values for the 
densities below) is more similar to those found for other Hii 
regions. 

Maps for both ratios are shown in Figure |6] According to 
the sulphur line ratio - detected over the whole field - densities 
range from very low values, of the order of the low density limit, 
in a region of about 5" x 5" in the upper right corner of the 
field just above the Complex #3, to 790 cm at the peak of 
the emission in the cluster associated with the H ii region, with 
a mean (median) over the field of view of ~ 130(90) cm 
The rest of the H ii regions still present high densities (of about 
400 cm as a whole, 480 cm in the Hii-2, (i.e. the up- 
per area of the giant Hii region, Kobulnickv et al. 1997). This 
agrees well with the value estim ated from long-slit measure- 
ments (ILopez-Sanchez et"aL 1 12007,) . The tail and the region as- 
sociated with the cluster UV-1 present intermediate values (of 
about 200 cm "^). 

The argon line ratio is used to sample the densest regions. 
The map for this ratio was somewhat noisier and allowed an es- 
timation of the electron density only in the giant H ii region. The 
densities derived from this line ratio are comparatively higher, 
with a mean (median) of 3400 (3150) cm As it happened in 
the case of the extinction, the peak of electron density according 
to this line ratio is offset by ~ 0'.'7 - 1'.' 1 towards the north-west 
with respect of the peak of continuum emission at Complex #1. 

We created higher S/N ratio spectra by co-addition of 3 x 
3 spaxel apertures associated with certain characteristic regions 
(i.e. the clusters at the core of the systems, C1+C2, and the re- 
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gions Hn-2, Hii-1 and UV-1 of iKobulnickv et alJ ([1993)). The 
largest values are measured around the core (i.e. C1+C2) where 
the [Ar IV] electron density can be as high as 6200 cm As 
we move further away from this region, the measured electron 
density becomes lower. Thus, H ii-2 presents similar, although 
slightly lower, densities (~ 6100 cm "^), followed by Hii-1 with 
~ 4200 cm and UVl with ~ 3300 cm These values agree, 
within the errors, with those reported in Lopez-Sanc hez et al.. 
( l2007h for similar apertures. 

An interesting point arises when the different density val- 
ues derived for the integrated spectrum and for each individual 
spaxel/aperture are compared (180 cm"-', and up to 790 cm"-', re- 
spectively when using the sulphur line ratio). The covered f.o.v. 
(~210pcxl35 pc) is comparable to the linear scales that one can 
resolve from the ground at distances of ~40 Mpc (or z ~ 0.01). 
Such a comparison illustrates how aperture effects can cause im- 
portant underestimation of the electron density in the H ii regions 
in starbursts at such distances, or further away. 



3.4. Ionization structure, excitation sources and nitrogen 
enhancement 

The ionization structure of the interstellar medium can be stud- 
ied by means of diagnostic diagrams. Different areas of a given 
diagram are explained by different ionization mechanisms. In 
the optical spectra l ran ge, the m ost wid ely used are probably 
those proposed by [ Bald win et alj (|1981|) and later reviewed by 
IVeilleux & Osterbrijckl d 19870 . the so-called BPT diagrams. In 
Figure [2[ the maps for the three available line ratios involved in 
these diagrams - namely [Nii]^6584/Hq', [S ii]^^67 17,673 1/Ha, 
[Oiii]/l5007/HjS - are shown on a logarithmic scale. This figure 
shows that the ionization structure in the central region of this 
galaxy is complex. Not only do the line ratios not show a uni- 
form distribution, but the structure changes depending on the 
particular line ratio. 

Both the [Oiii]^5007/H/3 and the [S ii]/1/16717,6731/Hq' line 
ratios display a gradient away from the peak of emission at 
Complex #1 and with a structure that follows that of the ion- 
ized gas. Thus the [S ii]^^67 17,673 l/Ha ([Oiii]/15007/Hj6) ra- 
tio is smallest (largest) at Complex #1, presents somewhat in- 
termediate values in the two tongue-shaped extensions and the 
Complex #2 and is relatively high (low) in the rest of the field, 
with a secondary minimum (maximum) at ~ [-3'.'5, -2'.'5], the 
position of a secondary peak in the Ha emission. This is coinci- 
dent with the structure presented in dCalzetti et alJl2004[) . 

The [N n]/}6584/HQ' line ratio, however, display a differ- 
ent structure. While in the right half of the FLAMES field 
of view, the behavior is quite similar to the one observed for 
the [S ii]/l/l6717,6731/Ha' ratio (i.e. values relatively high, lo- 
cal minimum at ~ [-3'.'5, -2'.'5]), the left half, dominated by 
the giant H ii region, displays a completely different pattern. The 
lowest values are associated with Complex #2 and the southern 
extension, rather than Complex #1, and the [Nii]/16584/Ha line 
ratios are highest at ~ [4'.'0, 2'.'0]. 

In Figure [8| we show the position of each spaxel of the 
FLAMES field of view, as well as for the integrated spectrum, 
after co-adding all the spaxels in the BPT diagnostic diagrams 
together with the borders that separate Hii region-like ioniza- 
tion from ionization by other mechanism s according to sev- 
eral authors ( Veilleux & Osterbrock 1987 : iKewlev et alj[200ll 
iKauffmannet al.1 i2003l: iStasinska et all 12006"). We also show 
the prediction s for models of photo-ionization caused by stars 
dPopita e t aL 2006) that take into account the effect of the stellar 
winds on the dynamical evolution of the region. In these mod- 
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Fig. 7. Emission line ratio maps. Up: log ([Nn]/}6584/Ha. 
Middle: log ([S n]/l/i67 17,673 1/Ha). Bottom: log 
([Oiii]/I5007/Fly6). The position of the three main peaks of 
continuum emission are shown for reference. 



els, the ionization parameter is replaced by a new variable R that 
depends on the mass of the ionizing cluster and the pressure of 
the interstellar medium {R - (Mc//M0)/(Po/k), with Po/k mea- 
sured in cm^ K). Also, the predictions for shocks models for a 
LMC metallicity are included. Given the relatively low metal- 
licity of NGC 5253, these are the most appropriate ones. They 
were calculated assuming a A^e = 1 cm"^^ and cover an ample 
range of magnetic parameters, B, and shock velocities, Vj, (see 
Allen et al. 2008, for details). 

As demonstrated for the electron density, these diagrams il- 
lustrate very clearly how resolution effects can influence the 
measured line ratios. Values derived for individual spaxels 
cover a range of ~0.5, 1.0 and 0.5 dex for the [Nii]/16584/Ha', 
[Sii]/1/16717,6731/Hq', and [Oiii]^5007/H/3 line ratios respec- 
tively, with mean values similar to the integrated values (-1.15, 
0.92, and 0.74). This is particularly relevant when interpreting 
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Fig. 8. Position of the individua l spaxels in NGC 5253 in th e 
diagnostic diagrams proposed bv lVeilleux & Osterbrockl(ll987h . 
Data points above and below the 3-cr level of the first-degree 
polynomial fit in the [Nii].16584/Hq' vs. [S ii]^^6717,6731/Ha' 
diagram (see text for details) have been represented with blue di- 
amonds and yellow crosses respectively. The solid curves show 
the empirical borders found by Veilleux & Osterbrock between 
ionization caused by different mechanisms, whil e the dotted 
lines show the theoretical borders proposed by Kewlev et al.l 
(l200l to delimit the area where the line ratios can be explained 
by star formation. Black dashed and dot-dashed lines s how the 
revised borders by KaufFmann et al. (2003) and Stasirisk a et alJ 
( |2006|) . respectively. Solid horizontal lines show the border be- 
tween classical Seyfert galaxies (above the line) and LINERs 
(below the line), at [Oiii]/t5007/HyS=3. Long-dashed lines in or- 
ange represent the predictions from models of photo-ionization 
by stars presented in Dopita et al. (2006) for R = and ages 
of 2 - 6 Myr The g rids represent the p redictions for the shock 
models presented in lAllen et aL with and without pre- 

cursor, assuming the abundances of the LMC and Ne = I cm"^. 
Lines of constant magnetic parameter (blue) cover a range of 
B = 0.0001 - 10 yuG, while lines of constant shock velocity (red) 
cover a range of Vj = 150 - 500 km s"'. Green asterisks mark 
the derived values for the integrated spectrum of the whole IFU 
field. 



the ionization mechanisms in galaxies at larger distances where 
the spectrum can sample a region with a range in ionization prop- 
erties. This loss of spatial resolution thus 'smears' the determi- 
nation of the ionization mechanism by a set of line ratios. Even 
if this given set of line ratios is typical of photoionization caused 
by stars, it is not possible to exclude some contribution due to 
other mechanisms at scales unresolved by the particular obser- 
vations. 



Regarding the individual measurements, although all line ra- 
tios are within the typical values expected for an H ii region-like 
ionization, two differences between these diagrams arise. The 
first one is that the diagram involving the [S ii]/l/l67 17,673 1/Ha 
line ratio indicates a somewhat higher ionization degree than the 
one involving the [Nii]/16584/Hq' line ratio. That is: values for 
the diagram involving the [S ii]/l/l6717,6731/Ha line ratio are 
at the limit of what can be expla ined by pure photo- ionization 
in an H ii region according to the iKewlev et alj ( 1200 Ih theoreti- 
cal borders. On the contrary, most of the data points in the di- 
agram involving the [Nii]/16584/Ha are clearly in the area as- 
sociated to photoionization caused by stars. A comparison with 
the predictions of the models for metallicities similar to the one 
of NGC 5253 shows how the measured line ratios present in- 
termediate values between those predicted by ionization caused 
by shocks and those by pure stellar photoionization. That this is 
exactly what one would expect if shocks caused by the mechan- 
ical input from stellar winds or supernovae within the starburst 
were contributing to the observed spectra. Also, this comparison 
supports previous studies that show how models of photoioniza- 
tion caused by stars underpredict the [O Iii1/150 07/H8 line ratios , 
specially in the lo w-metallicity cases (IBrinchmann et al.1 120081 : 
iDopita et al.ll2006l) . 
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Fig. 9. [Sn]/l/16717,6731/HQ' vs. [Nii]^6584/Ha' diagram. The 
red box includes the data points utilized to determine the 
first-degree polynomial fit (dark blue continuous line) between 
the [Sn]^/1671 7,673 1/HcK and [Nii]^6584/Hq' line ratios. Data 
points above and below the 3-cr level (light blue continuous 
line) have been represented with blue diamonds and yellow 
crosses respectively. The expected line ratios for H n regions for 
Z = O.2Z0 and Z - QAZq with ages in the r ange of 2-6 Myr and 
R = Q according to the models presented in Dopita et al.l (l2006h 
are over-plotted with dashed oranges lines for comparison. 



The second diference is the distribution of the data 
points in these diagrams. While the data points in the 
[Sn]/l/i67 17,673 l/Hff vs. [Oiii]/I5007/Hje diagram form a se- 
quence, data points in the [N ii]/}6584/Ha vs. [Oiii]/I5007/Hj6 
diagram are distributed in two groups: a sequence similar to 
the one in the [S ii]^^6717,6731/Ha' vs. [Oiii]^5007/Hy3 dia- 
gram and a cloud of data points above that sequence with larger 
[Nii]/16584/Hq'. This result can be interpreted either by local 
variations in the relative abundances or by changes in the ion- 
ization parameter. Here we will explore the first option, which 
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Fig. 10. Spaxels with [Nii]/16584/Hq' line ratio above the 3-cr 
level of the fit presented in Figure |9] on top of our He emission 
line map. The size of the white circles is proportional to the ex- 
cess in the [Nii]/16584/Hq' line ratio. Dashed and dotted lines 
delimit the apertures utilized to extract the spectra with Wolf- 
Rayet and nebular Hen features displayed in Figures [TT] and [12] 
(see sections [33] and [33] l. The position of the three peaks of con- 
tinuum emission in the map are shown by crosses for reference. 



is the most accepted explanation (e.g. iLopez-Sanchez et"al] 
l2007l and references therein) and is supported by the rela- 
tively constant ionization paramete r found in specific areas 
via long-slit (log([/) —3, Kobulnic kv et all [1997). Long-slit 
measurements in specific areas of this field have shown out 
how this galaxy present some region s with an over-abundanc e 
of nitrogen (e.g. IWalsh & Rovl 119891; iKobulnickv etal]|1997l) . 
For our measured line ratios and using expression (22) 
in [Perez-Montero & Contini (2009), we measure a range in 
\og{NIO) of -0.70 to -1.46. Here we will assume that this 
over-abundance is the cause of our excess in the [Nn]/i6584/HQ' 
line ratio and will use this excess to precisely delimit the area 
presenting this over-abundance. To this aim, we placed the in- 
formation of each of the spaxels in the [S ii]/1/16717,6731/Hq' 
vs. [Nn]/16584/Ha' diagram, which better separates the two 
different groups described above. This is presented in Figure 
|9] We have assumed that in the so-called un-polluted ar- 
eas, the [Nii]^6584/Hq' and [S ii]^^6717,6731/Hq' follow a 
linear relation. This is a reasonable assumption since the 
[Nn]/l6584/[Sii]/l/l6717,6731 line ratio has a low dependence 
with the abundance and the properties of the ionizing radiation 
field dKewlev & Dopita 2002). This standard relation was deter- 
mined by fitting a first-degree polynomial to the data points with 
[Sn]/l^67 17,673 1/Ha'> -0.8 (indicated in Figure |9] with a red 
box). Those spaxels whose [Nii]/16584/Hq' line ratio was in ex- 
cess of more than 3-cr from the relation determined by this fit, 
have been identified as having an [N IIl/Haexcess, and are iden- 
tified by diamonds in Figure |9] As can be seen from this figure, 
there are a number of spaxels where this excess is much above 
the standard relation. 

The data points thus identified with [Nii]/16584/Hq' excess 
are shown as a map in Figure [10] where the location and magni- 
tude of the [Nii]/16584/Hff excess is indicated by white circles, 
whose size is proportional to the size of the [Nii]/16584/Hq' ex- 
cess. This figure can be interpreted as a snapshot in the pollution 
process of the interstellar medium by the SSCs in the central area 
of NGC 5253. The pollution is affecting almost the whole giant 
H II region. The largest values are found at ~ r.'5 towards the 
north-west of the peak at Complex #1 . Then, the quantity of extra 



nitrogen decreases outwards following the two tongue-shaped 
extensions towards the north-west and south-east. This is con- 
sistent with the HST observations of Kobulnickv et al. (199^ 
who found nitrogen enrichment in their Hii-1 and Hii-2, while 
the N/0 ratio in UV-1 (see Table [3]) was typical for metal -poor 
galaxies. 

3.5. Wolf-Rayet features 

Wolf-Rayet (W-R) stars are very bright objects with strong broad 
emission lines in their spectra. They are classified as WN (those 
with strong lines of helium and nitrogen) and WC (those with 
strong lines of helium, carbon and oxygen) and are understood 
as the result of the evolution of massive O stars. As they evolve, 
they loose a significant amount of their mass via stellar winds 
showing the products of the CNO-burning first - identified as 
WN stars - and the He-burning afterwards - as WC stars ( Cor^ 
|1976|) . The presence of W-R stars can be recognized via the W- 
R bumps around /14650 A (i.e. the blue bump, characteristic of 
WN stars) and /15808 A (i.e. the red bump, characteristic of WC 
stars, but not covered by the present data). 

ISchaerer et al.l (1 19971 Il999h carried out a thorough search 
and characterization of the W-R population in NGC 5253. They 
detected W-R features (both WN and WC) at the peak of emis- 
sion in the optical (our Complex #1) and the ultraviolet (our 
Complex #2). Due to the spatial coinc idence of these dete ctions 
and the N-enriched regions found bv IWalsh & Rovl (1198 9) and 
Kobulnickv et al. ( 1997) - at least in the case of the nucleus - 
they suggested that these W-R stars could be the cause of this en- 
hancement. This is supported observationally by similar findings 
in other W-R galaxies. For example, a recent survey using Sloan 
data of W-R galaxies in the low-redshift Universe has shown that 
galaxies belonging to this group present an ele vated N/O ratio, in 
comparison with similar non-W-R galaxies jBrinchmann et al.1 
2008). Other suggested possibilities to cause the enrichment in 
nitrogen include planetary nebulae, O star wi nds. He-deficient 
W-R star winds, and luminous blue variables (IKobulnickv et al.l 
1997.). 

Here, we characterize the W-R population in NGC 5253 and 
explore the hypothesis of W-R stars as the cause of the nitro- 
gen enhancement by using the 2D spectral information provided 
by the present data. In the previous section we have delimited 
very precisely the area that presents nitrogen enhancement. In a 
same manner, it is possible to look for and localize the areas that 
present W-R emission. Note that due to the continuous sampling 
of the present data this can be done in a completely unbiased 
way. 

We visually inspected each spectrum looking for the more 
prominent W-R features in the blue bump (i.e. Niii/14640 and 
He 11/14686). The areas where these features have been found 
are marked in Figure [TO] with dashed lines. The co-added and 
extracted spectra of each individual area appear in Figure [TTIto- 
gether with a reference spectrum, free of W-R features, made by 
co-adding 20 spaxels in the upper right corner of the FLAMES 
field of view. 

We confirm the detection of W-R features associated with 
the nucleus and the brightest cluster in the ultraviolet, UV-1 (our 
W-R 3). In the same manner, we also detect a broad He ii line as- 
sociated with the north-west and south-east extensions (i.e. Hii- 
2 and W-R 2, respectively). In addition, there are three more 
areas which present W-R features, called W-R 1, W-R 4 and W- 
R 5, relatively far (~ 58 - 83 pc) from the main area of activity. 
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Fig. 11. Spectra showing Wolf-Rayet features. The first spectrum has been extracted from a low surface brightness area in the upper 
right corner of the FLAMES field of view and is presented here as reference. The positions of the nebular emission lines have been 
indicated with blue ticks and labels, while those corresponding to Wolf-Rayet features appear in red. The position and extent of the 
regions are indicated on Figure [TOl 



Interestingly, two of these regions (W-R 4 and W-R 5) present a 
narrow nebular He ii on top of the broad W-R feature. 

The short phase of W-R stars during star evolution makes 
their detection a very precise method fo r estimating the age o f 
a given stellar population. According to iLeitherer et al.l (Il999h . 
typically an instantaneous starburst shows these features at ages 
of ~ 3 - 6 Myr for metallicites of Z = 0.004 - 0.008, similar to 
the one in NGC 5253. Thus, very young stellar clusters must 
be associated with the areas that display these W-R features. 
We compared the p ositions of ou r det ections with the catalogue 
of clusters given by iHarris et"an (1200 4) and compiled in Table 
m All our regions, except W-R 5, are associated with one (or 
several) young (i.e < 5 Myr) star cluster(s). Regarding W-R 5, 
[Harris et aH (l2004l) do not report any cluster associated with that 
area. 

Cluster ages in Har ris et al.l (j2004) were estimated using both 
broadband photometry and the Ha equivalent width. We esti- 



mated the ages by means of two indicators: the ratio between 
the number of W-R and O stars; and the H/3 equivalent width. 
The ratio between the number of W-R and O stars was estimated 
from F(bb)/F(Hje), where F(bb) and F(Hj6) are the flux in the blue 
bump (measured with splot) and in H/? respectively and using 
the relation proposed bv Schlegel et al. (1998). Uncertainties are 
large, mainly due to the difficulty to define the continuum and 
to avoid the contamination of the nebular emission lines when 
measuring F(bb) but indicates a range in log(WR/(WR+0)) of 
2.0 to -1.2 (see TableH] column 4). There is a higher pro- 
portion of W-R stars in W-R 5, W-R 1 and W-R 4 and somewhat 
lower in those areas associated with the giant H ii region. The 
HyS equivalent widths are extremely high, consistent again with 
the expected youth of the stellar population. Predicted ages from 
these two age tracers are reported in Table |4] co lumns 5 and 7. 
They were estimated by using STARBURST99 dLeitherer et al.l 
Il999h . assuming an instantaneous burst of Z = 0.008, an up- 
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Harris et al 


(2004) 




This work 






Region 


Clusters 


Age 


log(WR/(WR+0)) 


AgewRyo 


EW(H/3) 


Age^a 




Number 


(Mwr) 




(Mvr) 


(A) 


(Mvr) 


Nucleus 


1 


3 


-1.99 


2.9 


245 


2.7 


Hn-2 


20 


3 


-1.93 


2.9 


320 


2.4 


W-R 1 


13 


4 


-1.34 


3.2/4.8 


131 


3.4 


W-R 2 


23 


5 


-1.64 


3.0/5.5 


180 


3.2 


W-R 3 


4,24,25 


1-5 


-1.55 


3.0/5.3 


128 


3.4 


W-R 4 


6,21 


4 


-1.35 


3.2/4.8 


131 


3.4 


W-R 5 






-1.17 


3.4 


86 


4.7 



per mass limit of M„p=100 and a Salpeter-type Initial Mass 
Function. The two age tracers give co nsistent age pred ictions 
and in agreement with those reported in [Harris et alj d200 4). 

The distribution of the W-R features, in an area of about 
100 pcxlOO pc, much larger than the one polluted with nitro- 
gen, suggests that all the detected W-R stars are not, in general, 
the cause of this pollution. Since the N-enrichment appears to 
be associated with the pair of clusters in the core and, given the 
position of maximum, most probably with the obscured SSC C2, 
the best W-R star candidates to be the cause of this enrichment 
are those corresponding to our Nuc aperture, and perhaps also 
the H II-2 and W-R 2 regions. 

3.6. Nebular He ii and helium abundance 

The hypothesis that the W-R population is the cause of the ni- 
trogen enrichment in NGC 52 53 requires an enhancement of 
the helium a bundance too (e.g . Schaererlll996 ). This is nicely 
illustrated in Kob ulnickv et all (1 1997 ) where different linear re- 
lations between the nitrogen and helium abundances (N/H and 
He/H) are presented according to different scenarios of nitrogen 
enrichment (W-Rs, PNe, etc.). The only scenario able to explain 
an extra quantity of nitrogen in the ISM without any extra he- 
lium counterpart would be the one where this nitrogen is caused 
during the late O-star wind phase. 

As in previous sections, we can measure at each spaxel 
the total helium abundance and compare it with that for nitro- 
gen. Since lines like [Oii]/l/l3726,3728 did not fall in the cov- 
ered spectral range, we did not determine the absolute nitrogen 
abundance. Instead, w e used t he me an of the abundances de- 
termined by Kobulnic kv et alJ (Il997h for their H ii- 1 and H ii- 
2 (N/H~ 2.0 X 10"^) to estimate how much helium would be 
needed in the enriched areas, if the extra nitrogen were caused 
by W-Rs (i.e. He/H~ 0.12). For the non N-enriched areas, we 
can use the measurement at UV-1 (N/H~ 0.7 x 10"^) which re- 
quires He/H~ 0.09. Helium abundance can be determined as: 



log(Hei:^6678/Ha) 



He/H = icf X (He^/H^ + He^^/H^) 



(1) 



where icf is a correction factor due to the presence of neu- 
tral helium. We assumed icf ~ 1.0, which is consistent with 
the predictions of photoio nization models for our m easured 
[Om]i5007/Hy6 line ratios (iHolovatvv & Melekhl 12002'). Since 
the He 1/16678 was detected in every spaxel of the FLAMES field 
of view and with good S/N, for the purpose of this work, we de- 
termined - He^/H^ from the He i/16678/HQf line ratio using 
the expression 
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Fig. 13. log(He I/16678/Ha) line ratio map. The position of the 
three peaks in the continuum map are shown for reference. 



y 



+ = 2.58f°-2^(2.87Hei.l6678/Ha) 



(2) 



where t is the electron temperature in units of 10"* K ('Pag el et akl 
1992). As in section [331 we assumed T,, - 11 650 K. Figure [13] 
shows the 2D structure of the He i/l6678/Ha line ratio. It is rel- 
atively uniform with the exception of some spaxels in the upper 
right corner, close to Complex #3. This area, relatively far from 
the main photo-ionization source and with low [Om]/l5007/H/3 
line ratio would be the only region where one can expect a sub- 
stantial contribution of neutral helium. 

Figure [14] presents the derived He^ abundances vs. the 
[Oiii]/I5007/Hy6 line ratio. With the exception of the data cor- 
responding to the right upper corner of the FLAMES f.o.v., 
He^/H^ range between 0.075 and 0.090, being higher in the 
higher excitation zones (i.e. the giant Hii region). These val- 
ues are in agreem ent with previous measurements of He^/H^ 
in specific a reas jPagel et al.[ [1992[ : [Kobulnickv et alj [1997[ : 
Walsh & Rov 1989). They are consistent with a scenario with- 
out extra N-enrichment and still far, by a factor ~ 1 .3 - 1 .7, from 
the required ~0.12 total helium abundance in the W-R scenario, 
in particular in the areas enriched with nitrogen. 

What about the He^^/H^, whose abundance can be deter- 
mined via the nebular He ii/l4686 emission line? This line turned 
out to be rather elusive. Campbell et al.. (.1986.) mentioned a pos- 
sible detection in their regions B and C. This result, however, has 
not been confirmed afterwards (see 'Lopez -Sanchez et al.[[2007l 
and references therein). As with the W-R features, we looked 
for the He n nebular line in each individual spectrum. Those that 
presented a spatial continuity were taken to define an area and 
were co-added before extracting. The selected areas are marked 
in Figure[TO]with dotted lines. The co-added and extracted spec- 
tra of each individual region appear in Figure [12] In addition to 
these regions, as mentioned in section 13.51 nebular He ii in the 
W-R 4 and W-R 5 has also been detected. 
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Fig. 12. Spectra showing nebular Hen, but no W-R features. The positions of the nebular emission lines have been indicated with 
blue ticks and labels, while those corresponding to Wolf-Rayet features appear in red. 
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Fig. 14. He"^ abundance vs. [Om]/I5007/H/? for each individual 
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Fig. 15. He++ abundance vs. [Oiii]i5007/H6 for each of the 
extracted spectra with Hen/14686 detection (dotted regions in 
Figure [Toll. 



These detections are, in general, neither associated with the 
area of nitrogen enhancement nor with those presenting W-Rs 
features. This lack of coincidence seems difficult to reconcile 
with a scenario where this enhancement, and the existence of 
He^^, share a common origin. Moreover, for the purpose of this 
work, we estimated the He ii abundances in this areas using: 



y 



0.084f"''*(He II/14686/H/3) 



(3) 



from'Pa gel etaP (Il992h . Derived values for the individual re- 
gions are shown Figure [15] as a function of the [Oiii]/I5007/Hj6 
line ratio. They range between 0.0001 and 0.0005. Although 
uncertainties are large, up to 0.0006, due to the weakness of 
the He 11/14686 emission line, these values are clearly very far 
from the values of ~ 0.030 - 0.050 required to bring the he- 
lium abundance up to ~0.12 in the W-R enrichment scenario. 
Given the depth and continuous mapping of the present data, 
we can clearly exclude the possibility of further detections of 
larger quantities of He^^ based on optical obse rvations. Thus the 
presen t data support the scenario suggested bv lKobuMckv et aP 
where the N-enrichment should arise during the late O- 
star wind phase. In view of the extra-nitrogen distribution and 
the extinction map, the only place where these larger quanti- 
ties of He^^ could be found (if they existed) is in Complex #1. 
However, they should be highly extinguished and would required 
a search for emission lines at longer wavelengths as for example 
Hen (7-10) at 21 891 A (iHora et alJll999l) . 

Since the nebular He n is not associated with the area show- 
ing N-enhancement, ther e is still the open question as to its ori- 
gin. iGarnett et alj (1199 Ih explored the different mechanisms ca- 
pable of producing this emission in extragalactic Hii regions. 
The first suggestion is photoionization by fast shocks. However, 
we have seen in section [34l that shocks do not appear to play a 
dominant role in the central parts of NGC 5253. Moreover, the 
measured log(He II/l4686/Hy6) are ~ -3.0 to -2.1, much lower 
than those predicted by shocks models with A^^ - Icm"^ and 
LMC abundances (~ -1.4 to -0.4, Allen et al. 200 8). 

Another possibility discussed bv lGarnett et al.l ( 119911) is hot 
(r<70000 K) stellar ionizing continua. This looks like a plau- 
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sible explanation for those cases where we had detected nebular 
He 11/14686 on top of the blue bump (WR 4 and WR 5). 

The last option would be photoi onization caused by X- rays. 
The only point sources detected by ISummers et alj (l2004h that 
fall in our f.o.v. are sources 17, 18, and 19. This last source ap- 
pears to be associated with the Complex #1 and thus, not related 
to this discussion (since no nebular He ii was detected for this 
region). Sources 18 and 17 could, however, be associated with 
Heii-l and Hen-4 detections, respectively. In particular, the lat- 
ter region coincides with the secondary peak of emission in the 
Ha image and is assoc iated with cluster 17 of the sample cata- 
logued by [Harris et alJ ('^04). No satisfactory explanation was 
found for the cause of the ionization at He ii-2 and He ii-3. 



3.7. Kinematics of tine ionized gas 

Slit observations in specific regions of this galaxy have demon- 
strated that the kinematics of the ionized gas is rather com- 
plex, with line profiles revealing asymmetric win gs (e.g. 
iMartin & Kennicut^ Il995t iLopez- Sanchez et alJ l2007h . These 
observations usually include the bright core of the galaxy. 
However, they might be biased since typically the long slits only 
sample specific regions selected by the particular slit placement. 
The present data permit a 2D spatially resolved analysis of the 
kinematics of the ionized gas in the central area of the galaxy 
to be performed, thus overcoming this drawback. We based our 
analysis on the strongest emission lines (i.e. mainly Ha, but also 
H/3, [Oiii].15007, [N ii].i6584, and [S ii]/l/l6717,6731) where the 
high S/N permits the line profiles to be fitted with a high degree 
of accuracy. 

In the following, we will present the results derived from 
Ha. Similar results were obtained from the Wfi and [Oiii]/I5007 
emission lines with differences of |Av|<2 km s"' in most cases 
and always between -5 and 5 km s"^ Results for the only emis- 
sion lines with remarkable differences in the velocity maps (i.e. 
[Nn].i6584 and [S ii].i.i6717,6731) will be presented, as well. 

Typical examples of the different profiles for the main emis- 
sion lines are shown in Figure [16] Lines are ordered by wave- 
length from bluer (lower) to red (upper) within each panel. The 
zero point in the abscissa axis corresponds to the measured sys- 
temic velocity which is defined as the average of the velocities 
derived from the main emission lines for the peak of the con- 
tinuum emission. We performed an independent fit for each of 
the brightest emission lines using MPFITEXPR (see section l231 
for details). In general, the line profiles of the individual spec- 
tra cannot be properly reproduced by a single component. A 
large percentage of them needed two (and even three) indepen- 
dent components to reproduce the observed profile reasonably 
well. We followed the approach of keeping the analysis as sim- 
ple as possible. Thus in those cases where both fits - the one 
with one component and the one with several components - re- 
produced equally well the line profile, we gave preference to 
the fit with one component. Examples of these fits are shown in 
Figure [T7] which contains the Ha emission line together with 
the total fit and the individual components overplot for the spax- 
els shown in Figure [161 Note how the fits for the spaxels in the 
central area of our f.o.v. present relatively larger residuals. These 
can be attributed to a low surface brightness broad extra compo- 
nent which would be the subject of a future work. Central wave- 
lengths were translated into heliocentric velocity taking into ac- 
count the radial velocity induced by the Earth's motion at the 
time of the observation which was evaluated using the IRAF 
task rvcorrect. Velocity dispersions were obtained from the 
measured FWHM after correcting for the instrumental width and 



thermal motions. The width of the thermal profile was derived 
assuming Tg^l 1 650 K which translates into a crther = yfkTJmu 
of ~11 km s ' for the hydrogen lines. The measured systemic 
velocity was 392 km s ' . This is slightly l ower than the one mea- 
sured from neutral hydrogen (407 km s~ MKoribalski et al.l2004l) 
according to NED. 

In Figure [18] we present the velocity fields for the three fit- 
ted components derived from the Ha emission line. We also in- 
cluded the velocity dispersion map for our broadest component. 
Corrected velocity dispersions for the two narrow components 
were, in general, subsonic and will not be shown here. The only 
exception would be an area at ~ [4'.'0, -2'.'0]. The line profiles 
in this area show how the narrow component presents a conti- 
nuity with the two narrow components at ~ [4'.'0, - I'.'O] as if it 
were the result of a strong blending of these two components. 
However, we were not able to properly deblend these two com- 
ponents by means of our line fitting technique. 

From Figure [18] it is clear that the movements of the ion- 
ized gas are far from simple rotation. For the discussion we will 
separate the emitting area into the zone corresponding to the gi- 
ant H 11 region and the rest. The zone of the giant H ii region, 
occupying roughly the left part of the FLAMES field of view, 
shows in the upper part, line profiles that can be explained by 
two components while in some spaxels of the lower part a third 
component was required. The area of the giant H ii region itself, 
which occupies an area of ~120 pcx60 pc, requires up to three 
components to properly reproduce the line profiles. They were 
named CI, C2 and C3, according to their relative fluxes. The 
first component (i.e. CI) accounts for the ~ 45 - 68% of the flux 
in Ha, depending on the considered spaxel. It is relatively nar- 
row and constant, with a Av ~10 km s"' over a distance of ~6" 
(~1 10 pc) with slightly bluer velocities in the spaxels associated 
with the edge of the upper and lower extensions. 

The second component (i.e. C2) accounts for the ~ 27-55% 
of the flux in Ha. It is symmetric with respect of an axis that goes 
through the Complex #1 in the north-south direction. In compari- 
son with the first component, it presents large velocity variations 
(i. e. Av ~ 70 km s"' over ~ 4'.'7 or ~ 86 pc) and is relatively 
broad (cr ~ 20 - 25 km s"' ). Low surface brightness broad com- 
ponents have been reported in starburst galaxies using a slit since 
more than a decade and have been the subject of several the- 
oretical ( e.g. Tenorio-Tagle et al.lll997t) a nd observatio nal (e.g. 
Castaned ret al.1 11990: Gonzalez-Delgado 'etai] |1994|) studies. 
They usually represent a small fraction (~3-20%) of the total 
Ha flux and have widths of cr ~700 km s Recently, 2D spec- 
troscopic analysis of very nearby starbursts have shown how lo- 
cally, the line width is somewhat smaller (cr ~50-170 km"', see 
Westmoquette et al. 2009, and references therein). This is un- 
ders tood in the context of the so-called Turbulent Mixing Layers 
(e.g. lSlavin et al ]fT993h . However, the high surface brightness of 
C2, together with the symmetry in the velocity field and its low 
widths made us to explore an alternative explanation for it (see 
section|4Ti). 



The third component (i.e. C3) is present in a small area of 
the field of about I'.'O diameter (i.e. 18 pc) and centered at about 
[4'.'0,-r.'0]. C3 has a width similar to the first component (CI) 
but displaced ~50 km s ' towards the blue. This third component 
appears in a location about r.'5 south-east to the position of the 
SSC complexes in the tongue-shaped extension described in sec- 
tion [TT] This area presents low values of extinction (see Figure 
|3]l and surface brightness (see Figure |2]l. This third component 
shows indications of more extended weaker emission which was 
not fitted. 
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We pointed out before that velocity maps derived from the 
[Nii]/16584 and [S ii]/l/l6717,6731 emission lines showed areas 
with important differences with respect to the one obtained from 
He This is the case for the giant Hii region. Figure [19] contains 
the measured velocity differences for the two brighter compo- 
nents (CI and C2) in the subset of the field corresponding to the 
giant H ii region while Figure |20] presents examples of the in- 
dependent fits for Ha, [Nii]/16584, and [S n]/}6717. Differences 
exist for both, the narrow (i.e. CI) and broad (i.e. C2), compo- 
nents in [Nii]/16584 and are relatively symmetric with respect 
to the peak of continuum emission Complex #1, but with oppo- 
site sign and slightly different directions (P.A.~ 50° and ~ 30° 
for CI and C2 respectively). Also, the range of velocity differ- 
ences, (vhq. - V[Nii]) is larger for C2 than for CI (~ 70 km s"' and 
~ 30 km s"', respectively). As illustrated in the right hand map 
of Figure [191 the broad component for the [S ii]/l/l67 17,6731 
emission lines also shows velocity differences with similar range 
(vHa - V[sii] ~ 60 km s"'), orientation and sign as in the case of 
the [N ii]/l6584 emission line. For the CI component, no signifi- 
cant differences were found. Measured differences in C3, with a 
mean and standard deviation of -6 + 9 km s"' and 5+4 km s"' 
for vuo - V[Nii] and vhq. - V[sii] respectively, do not appear to be 
significant. However, since C3 was only detected in seven spax- 
els (see Figure[T8] bottom left panel) this result has to be treated 
with caution. 

Si milar offsets has been dete cted in galactic H n regions like 
Orion (iGarcia-Diaz et al.l l2008h . but to our knowledge, this is 
the first time that maps with such offsets in velocity for dif- 
ferent emission lines in starbursts are presented. This can par- 
tially be caused by the fact that 2D-kinematic analysis of star- 
bursts, from dwarfs (e.g. G arcfa-Lore nzo et al. 2008)^to^ more 
extreme events like LlRGs (e.g. lAlonso-Herrero et al.lf2009l) . are 
usually based on fitting techniques that impose restrictions be- 
tween the Ha and [N n]/i6584 central wavelengths and thus, pre- 
venting from detecting such offsets. An example of work where 
the main emission lines are fitted independently is presented by 
IWestmoquette et al] (l2007h . However, since they only analyzed 
the kinematic for Ha, it is not possible to assess if they found 
different kinematics for the other emission lines. There are how- 
ever, some works that o ffer examples of offsets of t his kind us- 
ing a slit. In particular dLopez-Sanchez et alJl2007l) report also 
an offset between the [Nn]A6584 and the Ha emission line of 
~10 km s similar to what we have measured for CI in the 
upper part of our f.o.v. 

The second region of interest is located in the right (south- 
west) part of the FLAMES field of view. Emission there shows 
narrow lines with velocity dispersion dominated by the thermal 
width. In some areas (the north-east corner in Figure [T8l l. two 
narrow lines were needed to better reproduce the line profile. 
The primary component (i.e. CI) shows a symmetric velocity 
pattern with respect to the twin clusters associated with the peak 
of emission #3. A velocity gradient in the north-west to south- 
east direction is clear with a Av ~ 40 km over about 4'.'0 
(linear scale of ~75 pc). The secondary component traces a shell 
blue-shifted ~40 km s"' in the western corner (see Figure [T8l 
upper right panel). Note that C2, although relatively narrow, is 
a bit broader than the thermal width (Figure [18] bottom right 
panel). This component accounts for ~ 30 - 50% of the Ha flux 
in this area. No significant differences in the velocity fields and 
the velocity dispersion maps for the main emission lines have 
been found in this area. 



4. Discussion 

4.1. The giant H ii region 

The most interesting area of NGC 5253 in the present data cov- 
ers the left (north-east) part of the FLAMES field of view. In 
previous sections we have seen that this area is occupied by a 
giant H ii region which: i) harbors two very massive and young 
SSCs at its centre (i.e. Complex #1); ii) presents high levels of 
extinction, being larger in the upper part of the f.o.v.; iii) has 
high electron densities as traced by both the sulphur and the ar- 
gon line ratios, and again are also larger in the upper part of 
the f.o.v.; iv) presents an excess in the [Nn]/16584/Ha line ratio 
with respect to [S ii]/l/l67 17,673 1/Ha which, if interpreted as N- 
enrichment, indicates an outward gradient of extra nitrogen from 
a point at ~ 0'.'5 - I'.'O towards the north-west of the peak of 
continuum emission at Complex #1; v) presents W-R features, 
implying a young age for the harboured stellar population; vi) 
displays complex kinematics, not coincident in all the emission 
lines, that require a minimum of three components to reproduce 
the line profiles. How does all this evidence fit together into a 
coherent picture? 




Fig. 21. Sketch showing the different elements associated with 
the area of the giant Hii region. The two massive SSCs (stars 
in the circles) expel material (red and blue half moons) which 
encounter quiescent gas (grey shell). Regions of dust extinction 
are represented by the two wavy sheets. 

In Figure [21] we sketch a plausible scenario compatible with 
all these results. Here, the broad component (C2) would trace 
an outflow created by the two SSCs at Complex #1 while the 
two narrow components (CI and C3) would be caused by a shell 
of previously existing quiescent gas that has been reached by the 
ionization front. Different grades of grey in Figure|2T]in the shell 
represent the different densities observed in the upper and lower 
part of the FLAMES f.o.v., while two wavy sheets in two grades 
of grey have been used to represent the differences in extinction 
between these two halves. Due to this extinction distribution, in 
the upper (i.e. north-western) half of the H ii region, the observer 
cannot see the further part of the shell, while in the lower (i.e. 
south-eastern) half both parts are visible and are detected as a 
single broader component when approaching the vertex of the 
oval. 

As in section [373] we formed [S n]/16717/[Sii]/l6731 maps 
for the individual kinematic components and measured the mean 
and standard deviation for the sulphur line ratio in two 4 x 
4 spaxel squares sampling the upper and lower part of the 
giant Hii region. Although the standard deviations are large 
(~ 0.09), results support this sketch. While the broad compo- 
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nent presented similar line ratios in both areas (~ 1.07 im- 
plying densities of ~470 cm""^), the narrow component pre- 
sented somewhat lower line ratios in the upper part than in 
the lower one (~1.11 vs. ~1.25) implying densities for the 
shell of ~390 cm"^ and ~180 cm"^, respectively. Also, we 
created (noisier) [Nii]/16584/Ha', [S ii]^^6717,6731/Hq', and 
[Oiii]/I5007/Hj6 line maps and compared the relation between 
[Sii]^^67 17,673 1/Hq' and [Nii]^6584/Hq' for the individual 
components. The three fitted components present extra nitrogen 
in an area coinciding with the one derived from one-gaussian fit- 
ting. This differs from the findings for Mrk 996, a galaxy with 
several kinematically distinct components where only the broad 
one presented N-enrichment, with an abundance ~2Q times 
large r than the one for the narrow component (see iJames et"al] 
12009 ). Still, in NGC 5253, the N-enrichment in the broad com- 
ponent is larger than in the narrow one by a factor of ~ 1 .7 which 
is consistent with the scenario sketched above. Moreover, while 
the [Oiii]/I5007/Hj6 maps are relatively similar for both compo- 
nent^ (not shown), those associated to [N ii]/l6584/Ha' and spe- 
cially to (the more shock sensitive) [S n]/l/i67 17,673 l/Ho' line 
ratio display a relatively different ionization degree with larger 
values for CI than for C2. This is also consistent with the pre- 
sented scenario since a larger contribution due to shocks is ex- 
pected in the area where the outflowing material encounters the 
pre-existent gas. 

An interesting result of the previous section was the offsets 
derived for the velocities of the different species, in particular 
nitrogen and sulphur To our knowledge, this is the first time that 
maps showing this kind of offsets are reported in an starburst 
galaxy. Similar phenomena have already been reported in much 
closer regions of star formation. For example, observations in 
the Galactic Orion Nebula, a much less extreme event in terms 
of star formation, show how Ha and [O iii]/l5007 display similar 
velocities while [Nn1/I65 84 and [ S iiU/167 17,6731 are shifted 
by ~4-5 km s^' dGarcia-Diaz et al. 2008), an order of magnitude 
smaller than the shifts found for NGC 5253. Also, self-consistent 
dynamic models of steady i onization fronts poi nt towards the 
detection of such differences ("Hennev et alJ2005h . In the context 
of the scenario sketched in Figure [21] the offsets in C2 would fit 
if [Nii]/}6584 (and [S n]/l/i67 17,6731) traced the outer parts in 
an outflow which has a Hubble flow (i.e. velocity proportional to 
radius). 

Finally, an estimate of the time scales associated with the 
pollution process can be determined by using the velocity for 
the outflow derived in section [TTI Assuming that this traces the 
velocity of nitrogen contamination of the ISM, the detected pol- 
lution extending up to distances of ~60 pc took place over only 
~ 1.3 -1.7 Myr This supports the idea that the nitrogen dilution 
is a relatively fast process and is consistent with the shortage of 
observed systems presenting this kind of chemical inhomogene- 
ity. 

4.2. The area associated with the older steiiar dusters 

The rightmost (south-west) part of the FLAMES f.o.v. presents 
a different picture. We have seen that this region: i) is associ- 
ated with two relatively ol d (~70 and ~110 Myr) and massive 
(3 and 7x10"* Mq) clusters (iHarris et al.ll2004 : ii) presents mod- 
erate levels of extinction, being higher in the lower part of the 
FLAMES field of view; iii) has very low A^^ and lower than 
100 cm in the upper corner; iv) the Ha surface brightness is 



* C3 is not considered here, since the so-called map would be asso- 
ciated to < 7 spaxels, depending on the line ratio. 



Table 5. Integrated line ratios for the rightmost (south-west) part 
of the FLAMES f.o.v. 



Component 


log([Oiii]/H;e) 


log([Nn]/HQ') 


logCLSnl/Ha) 


Upper Blue 


0.69 


-1.20 


-0.85 


Upper Red 


0.52 


-1.09 


-0.67 


Lower 


0.54 


-1.09 


-0.75 


NGC 1952<^' 


0.92 


0.15 


0.18 



<"> From lOsterbrock & Ferlarig ( 120061) . 



very low (i.e. one and two orders of magnitude smaller than in 
the giant H ii region for the upper and lower portions respec- 
tively); v) displays two distinct kinematic components in the 
upper part of the f.o.v.. Noteworthy is that two supernov a rem- 
nant c andidates have been detected in the area (Labrie & Pritchetl 
l2006h . One of them (SOOl) appears very close in projection to 
the massive clusters at ~[-5'.'0,r.'0]. The second one (S002) is 
located, just outside of the FLAMES f.o.v., at the right upper 
corner. 




Fig. 22. Sketch showing the different elements associated with 
the right part of the FLAMES field of view. The two massive 
SSCs (stars in the circles) have already expelled their mate- 
rial and are seen as relatively quiescent ionized shells (shown 
in grey). Differences in extinction are represented by two wavy 
sheets. 

The first question to consider is if any of the kinematic 
components is related to the supernova remnant candidates. 
However, our measured velocity dispersions are much more 
lower than expansion velocities of typical supernova remnant s 
(i.e. NGC 1952, 1450 km s ' lOsterbrock & Ferlandl l2006h . 
Moreover, we created two integrated spectra for the upper 
and the lower part of the area. Line ratio for both compo- 
nents of the upper part and for the lower part were relatively 
similar (i.e. within ~0.1 dex, see Table |5]l and much lower 
than those expected for a su pernova remnant (e.g. NGC 1952, 
Osterbro ck & Ferlandll2006h . Thus, supernova remnants are not 
obviously the cause of the observed kinematics and physical 
properties of the ionized gas in this region. 

Instead, given that all the three components present similar 
line ratios - consistent with ionization caused by stars - and simi- 
lar line widths, a more plausible scenario would be that where all 
the components are part of a common picture. Given the age of 
the clusters, and the velocity differences between the three com- 
ponents, this area can be viewed as a snapshot of a more evolved 
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version of what is happening in the left part of the FLAMES 
f.o.v. The clusters have managed to clear out their environment. 
Only a broken shell made out of previously quiescent gas re- 
mains ionized by the remaining hot stars and moving away from 
the clusters with little evidence for high velocity outflow. Figure 
l22lpresents a sketch of the different elements associated with this 
area. 



5. Summary 

We present a thorough study of the ionized gas and its relation 
with the stellar population of NGC 5253 by mapping the central 
212 pcxl34 pc in a continuous and unbiased manner using with 
the ARGUS IFU unit of FLAMES. The analysis of the data have 
yield the following results. 

1. We obtained a 2D detailed map for the extinction suffered 
by the ionized gas, finding an offset of ~ 0'.'5 between the 
peak of the optical continuum and the extinction peak, in 
agreement with findings in the infrared. 

2. We compared the extinction suffered by gas and stars by 
defining ad hoc broad-band colours. We have shown the im- 
portance of using line-free filters when performing this com- 
parison and found that stars suffer less extinction than the 
ionized gas by a factor ~ 0.33, similar to the findings in other 
starburst galaxies. 

3. We derived A^^. sensitive line ratio maps. The one involv- 
ing the sulphur lines shows a gradient from 790 cm"^ at 
the peak of e mission in the giant Hii region described by 
ICalzetti et al.l (11997.) outwards. The argon line ratio is only 
detected in the area associated with the giant H ii region and 
traces the highest density (~ 4200 - 6200 cm"-*) regions. 

4. We studied the ionization structure by means of the maps of 
line ratios involved in the BPT diagrams. The spatial distri- 
bution of the [Sii]^^6717,6731/HQ'and [Oiii]/I5007/Hy3 line 
ratios follows that for the flux distribution of the ionized gas. 
On the contrary, the [N ii]/l6584/Ha' map shows a completely 
different structure. 

5. We evaluated the possible ionization mechanisms through 
the position of these line ratios in the diagnostic diagrams 
and comparing with the predictions of models. All our line 
ratios are compatible with photoionization caused by stars. 
The [S ii]/l/l67 17,673 1/Hq' indicated a somewhat higher ion- 
ization degree that might be evidence of some contribution of 
shocks to the measured line ratios. Part of the data in the di- 
agram involving the [Nn]/i6584/HQ' line ratio are distributed 
in a distinct cloud. This can be explained within the local 
N-enrichment scenario proposed for this galaxy. 

6. We delimited very precisely the area presenting local N- 
enrichment. It occupies the whole giant Hii, including the 
two extensions towards the upper and lower part of the 
FLAMES field of view, peaking at ~r.'5 from the peak of 
emission in the continuum and almost coincident (i.e. at 
~ 20 pc) with the peak of extinction. 

7. We located the areas that could contain Wolf-Rayet stars by 
looking for the blue bump. We confirmed the existence of W- 
R stars associated with the nucleus and the brightest cluster 
in the ultraviolet. W-R stars are distributed in a wider area 
than the one presenting N-enrichment and in a more irregular 
manner. We were able to identified one (or more) clusters 
with ages compatible with the existence of W-R stars in all 
but one (i.e. W-R 5) of our delineated regions with a W-R 
signature. 



8. If the scenario of N-enrichment caused by W-R stars turns 
out to be applicable, only the W-R detected at the core 
(Complex #1), and perhaps in the two extensions of the 
Giant H ii region, can be considered the cause of the local 
N-enrichment, according to the correlation of the spatial dis- 
tribution of W-R features and N-enrichment. 

9. We measured the He+ and He+^ abundances. He^/H^ is 
~0. 08-0. 09 in most of our field of view except for an area 
of ~2" X 3" in the upper right corner, far away from the 
main ionization source. We detected the nebular He ii/14686 
emission line in areas not coincident, in general, with those 
presenting W-R features, nor with the one presenting N- 
enrichment. Abundances in Hen were always <0.0005. 
Given the depth and unbiased mapping of the present data, 
we can exclude the possibility of further detections of larger 
quantities of He^^ based on optical observations in the nu- 
clear region of NGC 5253. This result is difficult to reconcile 
with the scenario of N-enrichment caused by W-R stars and 
favours a suggestion where the N-enrichment arises during 
the late O-star wind phase. 

10. We studied the kinematics of the ionized gas by using ve- 
locity fields and velocity dispersion maps for the main emis- 
sion lines. We needed up to three components to properly 
reproduce the line profiles. In particular, one of the compo- 
nents associated with the Giant Hii region presents super- 
sonic widths and [Nii]/16584 and [S n]/l/i67 17,6731 emis- 
sion lines shifted up to 40 km s ' with respect to Hot. Also, 
one of the narrow components shows velocity offsets in the 
[N n]/l6584 line of up to 20 km s"' . This is the first time that 
maps providing such offsets for a starburst galaxy have been 
presented. 

1 1 . We provide a scenario for the event occurring at the Giant 
Hii region. The two SSCs are producing an outflow that en- 
counters previously existing quiescent gas. The scenario is 
consistent with the measured extinction structure, electron 
densities and kinematics. 

12. We explain the different elements in the right (south-west) 
part of the FLAMES field of view as a more evolved stage 
of a similar scenario where the clusters have now cleared 
their local environment. This is supported by the low electron 
densities and Ha surface brightness as well as the kinematics 
in this area. 
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Fig. 16. Profiles of the main emission lines for representative spectra. The distribution within the FLAMES f.o.v. has been roughly 
retained. That is: spaxel (1,1) is located at the upper right comer of the FLAMES IFU and numbering increases towards the left 
and bottom. The middle row contains the profiles corresponding to the positions of the peaks of continuum emission (left: Complex 
#1; centre: Complex #2; right: Complex #3). Spectra in the left column correspond to spaxels associated with the giant Hn region 
(upper: extension towards the north-west; lower: extension towards the south-east), those in the central column to areas in the centre 
of the FLAMES f.o.v. and those in the right column to the area associated with Complex #3 (upper: position towards the north-west; 
lower: towards the south-east). 
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Fig. 17. Fits in Ha for the representative spectra shown in Figure [16] Flux is in arbitrary units. The observed spectra are presented 
in black. The first, second and third components have been displayed in blue, red, and yellow, respectively, while the underlying 
continuum and total fit are in violet and green. Residuals are shown as a dotted line below the profile. 
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Fig. 18. Kinematic information derived from the Ha emission line. Left: Velocity fields for the two narrow components. Right: 
Velocity field (upper) and velocity dispersion map (lower) for the broad component. 
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Fig. 19. Difference between the velocity fields derived from the Ha and the [N ii]/l6584 emission lines for the narrow (left) and broad 
(center) components as well as between those derived from the Ha and [S ii]/l/l67 17,6731 emission lines for the broad component 
(right). Only the section of the FLAMES field of view associated with the giant H ii region is shown. The "X" and the "O" symbols 
in the central panel show the position of the spaxels (17,4) and (17,9), repectively (see line fitting in Figure |20]|. 
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Fig. 20. Fits in Ha (left), [N ii]/l6584 (middle), and, [S ii]/l67 17 (right) for two representative spaxels showing offsets in velocity for 
the different emission lines. Flux is in arbitrary units. The observed spectra are presented in black. The first and second components 
have been displayed in blue and red, respectively, while the underlying continuum and total fit are in violet and green. Residuals are 
shown as a dotted line below the profile. 



